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We construct three families of general magnetostatic axisymmetric exact solutions of Einstein- 
Maxwell equations in spherical coordinates, prolate, and oblates. The solutions obtained are then 
presented in the system of generalized spheroidal coordinates which is a generalization of the previous 
systems. The method used to build such solutions is the well-known complex potential formalism 
proposed by Ernst, using as seed solutions vacuum solutions of the Einstein field equations. We show 
explicitly some particular solutions among them a magnetized Erez- Rosen solution and a magnetized 
Morgan-Morgan solution, which we interpret as the exterior gravitational field of a finite dislike 
source immersed in a magnetic field. From them we also construct using the well known "displace, 
cut and reflect" method exact solutions representing relativistic thin disks of infinite extension. 
We then analyze the motion of electrically charged test particles around these fields for equatorial 
circular orbits and we discuss their stability against radial perturbations. For magnetized Morgan- 
Morgan fields we find that inside of disk the presence of magnetic field provides the possibility of to 
find relativist charged particles moving in both prograde and retrograde direction. 
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I. INTRODUCTION 



Magnetic fields play an important role in the study of astrophysical objects such as neutron stars, white dwarfs, 
pulsars, black holes and galaxy formation. In fact, several observations show that there are various scenarios where 
the magnetic fields and general relativity can not be neglected. One of them is the presence of strong magnetic fields 
in active galactic nuclei [lH^] . These nuclei are known to produce more radiation than the rest of the entire galaxy and 
directly affect its structure and evolution. Another scenario is the production of relativistic collimated jets in the inner 
regions of accretion discs, which can be explained considering magneto-centrifugal mechanisms [El-Hol]. Also, magnetic 
fields are important in understanding the interplay between magnetic and thermal processes for strongly magnetic 
neutron stars pH4l3j . At least 10% of all neutron stars are born as magnetars, with magnetic fields above 10 14 G 
p^4l6j . Anal ytic al models that describe these astrophysical objects are often associated with solutions of Einsteins 
equations |17H22|. In the search for more realistic models for compact stellar systems, the energy-momentum tensor, 
the source of Einsteins equations, is modified by introducing more complex terms that take into account additional 
physical properties as, for example, electromagnetic fields [23|. 

Stationary or static axially symmetric exact solutions to the Einstein field equations representing relativistic thin 
disks are of great astrophysical importance since they can be used as models for certain galaxies, accretion disks, and 
the superposition of a black holes and a galaxy or an accretion disk as in the case of quasars. Static thin disks without 
radial pressure were first studied by Bonnor and Sackfield [88[ , and Morgan and Morgan [25[ , and with radial pressure 



by Morgan and Morgan [26j . Also thin disks with radial tension were considered [27[ • Several classes of exact solutions 
of the Einstein field equations corresponding to static thin disks with or without radial pressure have been obtained by 
different authors [28|, l30j-l3p . Rotating thin disks that can be considered as a source of a Kerr metric were presented 
by Bicak and Ledvinka [36| , while rotating disks with heat flow were were studied by Gonzalez and Letelier [37| . Disk 
sources for stationary axially symmetric spacetimes with electromagnetic fields, especially magnetic fields, are also of 
astrophysical importance in the study of neutron stars, white dwarfs and galaxy formation. In such situation one has 
to study the coupled Einstein-Maxwell equations. Thin disks in presence of electromagnetic field have been discussed 
as sources for Kerr-Newman fields [H, [33] , magnetostatic axisymmetric fields [40| , conformastationary metrics [4l[ , 
while models of electrovacuum static counterrotating dust disks were presented in (42j |. Charged perfect fluid disks 
were also studied by Vogt and Letelier (43L an d charged perfect fluid disks as sources of static and Taub-NUT-type 
spacetimes by Garcia- Reyes and Gonzalez [H, [45| . 

In all the above cases, the disks are obtained by an "inverse problem" approach, called by Synge the " g-method^ [jjjj]. 
The method works as follows: a solution of the vacuum Einstein equations is taken, such that there is a discontinuity 
in the derivatives of the metric tensor on the plane of the disk, and the energy-momentum tensor is obtained from the 
Einstein equations. The physical properties of the matter distribution are then studied by an analysis of the surface 
energy-momentum tensor so obtained. Another approach to generate disks is by solving the Einstein equations given 
a source (energy-momentum tensor). Essentially, they are obtained by solving a Riemann-Hilbert problem and are 
highly nontrivial [47H53l| . A review of this kind of disks solutions to the Einstein-Maxwell equations was presented by 
Klein in [EJ] . 

Motion of matter near compact stars and black holes has been discussed widely in literature. The interplay 
between gravitational and electromagnetic interaction is essential for characteristics of the motion, namely its stability 
properties. Motivation for these studies arises from the problem of motion and acceleration of matter (charged particles 
or dust grains) [55U58| . The study of the interaction between particles and electromagnetic fields in curved spacetimes 
is also of astrophysical interest, such is the case of strong synchrotron radiation emerging galactic cores, which can 
be explained admitting the existence in those regions of extended and very intense magnetic fields, interacting with 
ultrarelativistic electrons. Such magnetic fields could originate in the inner part of an accretion disc around the 
central black hole [5^, H3] ■ Also have been shown that the presence of a strong magnetic field provides the possibility 
of relativistic motion of counterrrotating matter 61 1. 

In the present work we construct three families of general magnetostatic axisymmetric exact solutions of Einstein- 
Maxwell equations in spherical coordinates, prolate, and oblates. The zero order of this solutions were investigated 
previously in the Ref. [42J, and although the models of thin disks constructed there satisfied all the energy conditions, 
the solutions are not asymptotically flat. Thus the solutions presented in this work are a generalization of solutions 
discussed in that reference, and, except the zero order, the new families of solutions are asymptotically flat but with 
total mass zero. Hence, all discussed solutions with magnetic field are either singular, not asymptotically flat or 
massless. Moreover, in all cases we build models of thin disk in which all the energy conditions are satisfied, which 
characterizes a matter distribution made of usual matter. In the case of oblates coordinates, the built solutions 



correspond to the magnetized version of the Morgan-Morgan vacuum solution 25j. These solutions are interpreted 
as the exterior gravitational field of a finite dislike matter distribution immersed in a magnetic field, and thus can be 
used to model disklike astrophysical objects such as galaxies, accretion disks, and certain stars in presence of magnetic 
fields. However, unlike the part material or the electrostatic case, we show that the source of the magnetic field is not 



3 



planar but of a different origin such as a remnants or fossil magnetic field [62], or can come from external sources, such 
as the presence of a nearby magnetars or neutron stars. On the other hand, even though realistic disklike sources have 
thickness, in first approximation these astrophysical objects can be considered to be very thin, e.g., in our Galaxy the 
radius of the disk is 10 kpc and its thickness is f kpc. In all cases we also construct relativistic models of infinite thin 
disks. We find for all the values of parameters that the surface energy density decreases rapidly which permits one can 
define a cut off radius and, in principle, to consider these disks as finite. In addition, we analyze the motion of charged 
test particles around these fields for circular equatorial orbits and we discuss the stability of this orbits against radial 
perturbations. For magnetized Morgan-Morgan fields we find that inside of disk the presence of magnetic field provides 
the possibility of to find relativist charged particles moving in both prograde and retrograde direction. Even though 
this result can be seen as merely theoretical, there are observational evidence of counterrotating matter components 
in certain types of galaxies |63l - l67| . 

The paper is organized as follows. In section II we make a review of the Ernst's method in the case of a magneto- 
static axisymmetric spacetime and we construct three families of solutions of Einstein-Maxwell equations in spherical 
coordinate, prolate and oblate, and then in generalized spheroidal coordinates which are a generalization of the pre- 
vious cases. We also show explicitly some particular solutions among them a magnetized Erez-Rosen solution and 
magnetized Morgan-Morgan solution, and we analyze limiting cases. 

In section III we present a summary of the procedure to obtain models of thin disks with a purely azimuthal 
pressure and currents for the Einstein-Maxwell equations. In particular, we obtain expressions for the surface energy- 
momentum tensor and the surface current density of the disks. In the case of the magnetized Morgan-Morgan solutions, 
these are interpreted as the exterior gravitational field of a finite dislike source immersed in a magnetic field. Then we 
construct using the well known "displace, cut and reflect" method exact solutions representing relativistic thin disks 
of infinite extension. 

In section IV we study for some particular solutions the equatorial circular motion of charged test particles and we 
also discuss their stability against radial perturbation using an extension of Rayleigh criteria of stability [68U70] | . In 
particular we analyze the circular velocity and the specific angular momentum of the particles. In order to compare 
the behavior of these physical quantities with other known magnetized solutions, we also study the electrogeodesic 
motion of test particles and their stability for a Kerr- type solution or magnetic dipole solution ;74] . Finally, in Section 
V we summarize and discuss the results obtained. 



II. GENERAL MAGNETIZED WEYL SOLUTIONS 

The simplest metric to describe a static axially symmetric spacetime is the Weyl's line element [73 



ds 2 = - e 2 ^dt 2 



[p 2 d<p 2 +e 2A (dp 2 +dz 2 )}, 



(1) 



where (t, tp, p, z) are the Weyl canonical coordinates, and ip and A are functions of the coordinates p and z only. For 
the coordinates we also use the notation (x , x , x 2 , x 3 ) = (t, ip, p, z). The vacuum Einstein-Maxwell equations, in 
geometrized units such that G = c = I, are given by 

R ab = 87rT Qb , (2a) 



Tab — ~r- 

4w 



F ab , b = 0, 



F ac F c — -g a bF c dF c 



(2b) 
(2c) 



F ab = A b ,a - A a , b , (2d) 

where T a b is the electromagnetic energy-momentum tensor, F a b is the electromagnetic field tensor, and A a = (</>, A, 0, 0) 
is the four-potential, where <fi is the electric potential and A the magnetic potential which are also functions of r and 
z only. Further, ( ) >a = d/dx a , and ( ) ;a means covariant derivate. 

For the metric ([T]) and in magnetostatic case, the Einstein-Maxwell equations are 



V • [p- 2 fVA] = 0, 



/V 2 / - V/ • V/ 

/.2 .1.2 



2p~'f A VA-VA, 
I 



A, p = p(iP 2 , p -ip 2 , z ) + -(A 2 , p 
2 

A, z = 2pip, p i/>, z +-A, P A, z f, 
9 



A 2 , z )f, 



(3a) 
(3b) 

(3c) 
(3d) 
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where / = e 2 ^ . The equations (|3a|) y (|3bl) are equivalent to [7l|, uM 

fA£ = CVS + 2**V$) • VS, (4a) 

/A$ = (V£ + 2$*V$) • V$, (4b) 

where A and V are the standard differential operators in cylindrical coordinates, / = e 2 ^, and £ and <£> are complex 

potentials which in the case static £ — £ * and for the case magnetostatic $* = — $ . The above equations are called 
Ernst equations. The metric functions are obtained via 

/ = £ + $$*, (5a) 



^(£ if +2M lC )(£ lC +24«; f )- ^ P 
where \/2( — p + iz, so that \/2d,Q = d. p — id. z , and the magnetic potential A is related to $ via 



A lC = -^(£ )f +2$*$ iC )(£, c +2$$* c )- ^$, c $; c , (5b) 



A c =^(I m <D) )f . (6) 



Taking £ as function of $, from (|4a[) and (j4b|) it follows that 

(Rrf + |$| 2 )^V$.V$ = 0, (7) 

and hence that £ is a lineal function of $. Using the boundary conditions £ — > 1 and $ — >• at infinity, we obtain 

£ = l-2g -1 $, (8) 
where g is a complex constant. With the change of variable 



then © implies 



*=(?TTy (10) 



and the Ernst equations read 

Kr-(i-99*)]v 2 e-2rve-ve (n) 

Then making £ = (1 — qq*) 1 ^ 2 S, the equation ([TT]) takes of form 

(|f-l)V 2 |=2fv|-Ve, (12) 

which is the Ernst equation in the vacuum [72| . So given a solution of the Einstein field equations in vacuum £ (seed 
solution) we can construct a solution of the Einstein-Maxwell field equations. A solutions for this equations is 

|=-e ia coth^, (13) 

where the function ip satisfies the Laplace's equation 

V 2 V> = 0. (14) 

The case a — corresponds to the well-known Weyl vacuum solutions and the magnetized solutions (taking $ 
imaginary) built from them are the magnetized Weyl solutions. The metric functions and magnetic potential are 
given by 

f = - -2, (15a) 

(1 + a) e-^ + (1 - a) 

A = A, (15b) 
A p = -bptp, z , (15c) 
A, z = b(4, P , (15d) 
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where a = y/1 + b 2 , being b is the parameter that controls the magnetic field, and A is the metric potential A 
corresponding to the seed solution, that is taking £ — e 2 ^ and $ = 0. 

In terms of generalized spheroidal coordinates (£,77) which are related to Weyl coordinates (p,z) by 

z = fcfty, 

with k and a constants, the equations for the magnetic potential (|15c[) - (|15d[) can be cast as 

A, t = -kb(l-rf)^ ini 
A, v = kb (e - a 2 ) fat . 

and for the function A 

1 -if 



A,, = 



A 



£ 2 - a 2 J] 2 



i (e - a 2 ) tl>* -£ (1 - V 2 ) j>* -2ij (e - a 2 ) ^ i>, v 
V if ~ ° 2 ) V3, 2 -r, (1 - if) fj 2 +2£ (1 - if) fa fa v 



(16a) 
(16b) 

(17a) 
(17b) 



(18a) 



(18b) 



When a = 1 we have the prolate coordinates, when a — i we have the oblate coordinates, and the case a — with 
i] = cos 9 and k = 1 corresponds to the spherical coordinates. As ip satisfies the Laplace's equation, the integrability 
of the systems (|17a|) - (Il8b|) is guaranteed. So we can choose any of these equations for find A or A. 



A. Solutions in spherical coordinates 

In terms of spherical coordinates (r,9) which are related to Weyl coordinates (p, z) by 

p = rsmd, z = rcos9, (19) 
with < r < 00 and — tt < 9 < it, the asymptotically flat general solution of Laplace's equation (|T4|) can be written 



Pn (COS t 
r n+l 



(20) 



n=0 



where c n are constants and P n (cos9) are the Legendre polynomials. The magnetic potential obtains from of any the 
equations (|17ap - (|17bl) in spherical coordinates 



A= A, g d9 = -b r 2 sm9ij;, r d9, (21) 

J IT J 7V 

where the integral limits are chosen by requiring that the function A to be regular on the axis of symmetry. So 

A = -b^c n {n + l)— I sm9P n (cos9)d9 
71=0 r ™ J* 

= &y^Cn(n+ 1)— / P n (cos9)d(cos9). 
With the change of variable y = cos 9 the above expression takes the form 

4 = ^c> + l)- / P n (y)dy 

n=0 T 



and using the identity [7£ 



P n (y)dy 



1 



2n+ 1 



(Pn+iiy) - P n -i(y)) , n>0, P_i=-1, 



(22) 
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we obtain 



I (2n + 1 r™ 



^ (2n+l) 

71=0 V 



Finally, using the identity 76 1 



2n + 1 

P n+ i(cos9) — P n ^i(cos9) = (cos 8P n (cos 9) — P n -i(cos8)) 

n + 1 

one finds that the magnetic potential can be written as 

oo 

A = b Yl ~Ti [cos 6LP n (cos 0) - P n _i(cos0)] . 

n=0 

The same above procedure is carried out to find A |73[ and we obtain 



(23) 



(24) 



(25) 



A = 



A, e d9 

ClCm {I + 1) (m + 1) 



E 

l,m=0 



(I + m + 2)r 



l+m+2 



[PP m - Pl+xP m +l] 



(26) 



Since r = y/ p 2 + z 2 the solutions have only a true singularity at the point p = 0, z — 0. Hence the solutions are 
regular on the axis of symmetry, except at the origin. Such as singularity could be interpreted as the inside some 
kind of black hole and the solutions as the gravitational field produced by such source. The Curzon-Scott-Szekeres 
black holes [77l [78j provides an example of this type of spacetimes. However, the problems associated with the 
complicated structure of the singularity in these space-times can be avoided by considering the solution to represent 
the exterior field of some finite source, e.g, representing the field of a disk (section III). In addition, since | cos(0)| < 1 
and jP^cos^)! < 1 [76| the magnetic potential A (|25t tends to zero at infinity only for n > 1. Hence, the solutions 
are asymptotically flat only for n > 1, i.e., if the coefficient cq vanishes. On the other hand, for the line element (p}, 
in terms of the spherical coordinates (|19|) . 



ds 2 



e W d t 2 + e -^[r 2 sin 2 0d<f + e 2A (dr 2 + r 2 d8 2 )] 



the total mass 176 



m 



— lim 

4"7T r— >-oo 



K a ' b dS ab 



(27) 



(28) 



(where S r is a coordinate sphere with radius r and K = dt is the static Killing vector) for an asymptotically flat 
solution is given by 



For a seed solution of the form ([2U)l we have 



777 



From (|15e 



Vv = i>,T 



lim (r 2 ip tr ). 



lim (r 2 7/>. r ) = — cq. 



(1 + a)e-^ ~ (1 - a)e^ 



(1 + a)e-^ + (1 - a)e^ 
so that using (|30|) and the fact that tp — > for r — > oo we obtain 

(1 + a)e-$ - (1 - a)e^ 



(29) 



(30) 



(31) 



lim (r tp^r) lim 



(1 + a)c-^ + (1 - a)e^ 



-ac . 



(32) 
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Now if one chooses cq — in order to guarantee asymptotical flatness this leads to m = for arbitrary coefficients 
Ci, C2, etc. These parameters contribute to higher multipole moments [80j but not to the mass. According to the 
positive mass theorem this imply that the spacetime must be singular or Minkowski space. So when the magnetic field 
is present, i.e., if b ^ 0, these solutions are massless in the asymptotically flat case. Moreover, since the expression 
for the mass (|28[) is invariant, this result is true for any coordinate system. Therefore, in presence of magnetic field 
these solutions are either singular, not asymptotically flat or massless. 

For n = we have that 

4 = ~, (33) 
r 

where we have chosen cq — m. Therefore, the exact solution of Einstein-Maxwell field equations for n = is 

= i ; t-^-, ; r, (34a) 

(f + a)e m ' r + (1 - a)e- m / r 

2 

771 

A = -— sin 2 6>, (34b) 
A = bm(cos6 + 1). (34c) 



This solution is the magnetized version of the Chazy-Curson vacuum solution [42, |81|, |82| ■ The solution is regular 
on the axis of symmetry but not is asymptotically flat. Therefore, this solution represents the exterior gravitational 
field of a source of finite mass in the presence of a magnetic field spread throughout the space (infinite) . 

We now consider the sum of the first and second terms n = and n = 2 of series (l20l) . In this case we have 



r 2r 

and the general solutions in spherical coordinates in this case is 



^ = ~ - £j(3 cos 2 9~l), (35) 



3 I 

A = bc (cos9 + 1) bc 2 ^ cos 6 sin 2 6, (36a) 

2 r 2 

A = -— sin 2 6 [3c 2 cos 2 9 (25c 2 cos 2 9 + 10c o r 2 - 14c 2 ) 



•sr 

+Aclr A - 6c Q c 2 r 2 + 3^] , 



(36b) 



and tp is given by ([T]). 



B. Solutions in prolate spheroidal coordinates 

In terms of the prolate coordinates (x, y) which are related to Weyl coordinates (p, z) by 

p 2 =fcV-I)(I-y 2 ), (37a) 
z = kxy, (37b) 

with x > 1, — I < y < 1, and k a constant, the asymptotically flat general solution of Laplace's equation (fT4"]l can be 
written as 

00 

$ = -^T,CnQ n (x)P n {y), (38) 
n=0 

where c„ are constants and the Q n (x) are the Legendre functions of the second kind 

y»W-2^ S !(2n + 2 S + I)! X ' ^ 

s— 

From (|17ap - f|17bl) in prolate coordinates, we obtain 

v 

A = -kb (x 2 - 1) CnQ' n {x) / P n (y)dy, (40) 
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where the integral limits are chosen by requiring that the function A to be regular on the axis of symmetry. Using 
the identities (|2"2"|) and (IM1) . we have that 



A = -kb (x 2 - 1) jT ^fjQn(x) [yPn(y) P n -i(v)]. 



(41) 



Since Q n {x) diverges for x — > 1, the resulting solutions are singular at x = 1, corresponding to p = 0, —1 < y < 1, 
i.e., on a part of the symmetric. In the vacuum these solutions have been interpreted by several authors as distorted 
static black holes [Hj]. Using the definition for Q n (x) (|39p . the recurrence relation 



(x 2 - l)Q' n (x) = (n + l)[Q„+i(x) - xQ„(x)] 



and the fact that Q n (oo) = we find that in the infinity (x — > oo) 



lim (x 2 - l)Q;(a;) 



-1, 71 = 

0, n > 1 



Therefore, the magnetic potential A (|4Tj) tends to zero at infinity, 

lim A = 0, 



(42) 



(43) 



(44) 



only for n > 1. Hence, the asymptotically flat solutions are those for which n > 1, i.e., if the coefficient Co vanishes. 
Therefore, as was discussed above, in presence of magnetic field these solutions are either singular, not asymptotically 
flat or massless. 
For n = we have 



%b = — In 

^ 2 



x- 1 
x + 1 



and the magnetic potential is 



Using 



A = -kbc (x 2 - 1) Q' (x) [yP (y) - P^y)} , 
- -kbc (z 2 -l) (y + l)Q' (x). 



h{x) = -In 



x + 1 
x- 1 

1 



(x 2 - 1)' 

we obtain 

A = kbc (y + 1). 
So the solution in prolate spheroidal coordinates for n — is 

^ _ 2(x 2 - l) c "/ 2 



(1 + a)(x + l) c ° + (1 - a)(x - l) c o ' 



A= ^ln 



1 



x 2 — y 2 



A = kbco(y + l). 



(45) 



(46) 



(47) 
(48) 

(49a) 

(49b) 
(49c) 



This solution is the magnetized version of the Zipoy- Voorhees vacuum solution [42|, [8J, |85| . For 6 = and Co = 1 
we have the Schwarzschild solution (8(| and the case b = and Co = 2 corresponds to the metric of Darmois [73| . 
The solution is regular on the axis of symmetry but is not asymptotically flat. Therefore, this solution also represents 
the gravitational field outside a source of finite mass in the presence of a magnetic field spread throughout the space 
(infinite). 

For the terms n — and n = 2 we have that 



V» = c i In 



x- 1 
x + 1 



\c2 (3y 2 - 1) 



-(3x 2 - l)ln 
4 





3 


(x + l) 









(50) 
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and for the second term of the summation n = 2 the magnetic potential is 



A(n = 2) = -kb (x 2 - 1) c 2 l -y (l - y 2 ) Q' 2 (x), 



= --kbc 2 y (l - y 2 ) 



3x (a; 2 -!) In 



x- 1 
x+1 



6x^-4 



Therefore the solution in prolate spheroidal coordinate in this case is 



A = fc6c (y + 1)- \kbc 2 y (l - y 2 ) 



3x (x 2 - 1) In (fL_l) + 6x 2 - 4 



A 



ci(^-l)(^-l)(9a:V-^ + l-^) 



In 



x- 1 
x + 1 



64 

+ 4c 2 x(2/ 2 - l)(27c 2 xV - 3c 2 a; 2 - 21c 2 ?, 2 + 5c 2 + 8c ) In ( 
16 \x + 1 

+ ^(c + c 2 ) 2 ln ^ ' ~ ' 



9 9 



"16 C2(2/2 ~ 1 )(" 12c2 y 2 + 27c 2zV - 3c 2 o; 2 + 4c 2 + 16c ), 



(51) 



(52) 



(53a) 



(53b) 



and ip is given by (|15a|) . In the absence of magnetic held 6=0 and Co = 1, we have the Erez-rosen vacuum metric 
[87j so that the case 6^0 corresponds to a magnetized Erez-rosen metric. 



C. Solutions in oblate spheroidal coordinates 

In terms of the oblate coordinates (u, v) which are related to Weyl coordinates (p, z) by 

p 2 =k 2 (u 2 + l)(l-v 2 ), 
z = kuv. 



(54a) 
(54b) 



with u > 0, —1 < v < 1, and k a constant, the asymptotically flat general solution of Laplace's equation (TT4"|) can be 
written as 



i> = - c n q n (u)P n (v), 



n=0 



where c n are constants and 

q„(w) = i n+1 Q n (iu). 

Again requiring that the solution to be regular on the axis of symmetry, it follows that 



V 

OO r. 

A = -kb (u 2 + l)Y, C X( U ) I P n (v)dv, 
n=0 i 



= -kb (u 2 + 2^T^(") [ P ™+^ ~ 

n=0 

and using the identities (|22|) y (|24|) we obtain the following expression for the magnetic potential 

OO 

A = -kb (u 2 + 1) ]T -xt^M - P B -i(v)], 

* — ' n + 1 

n=0 



(55) 



(56) 



(57) 



(58) 



and again i/j is given by (|15a|) . In the absence of magnetic held 6 = and for n even this solutions correspond 
to the Morgan-Morgan vacuum solutions [25| which represent the exterior gravitational field produced by a finite 
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disklike source. So we can call above solutions magnetized Morgan-Morgan solutions. Similarly to the prolate case, 
the asymptotically flat solutions are those for which n > 1, i.e., if the coefficient cq vanishes. Hence, as was discussed 
above, in presence of magnetic field these solutions are either singular, not asymptotically flat or massless. 
For n = we have 



and the potential magnetic is given by 



with 



1p = — Cq COt 1 (u), 

A = -kbc (u 2 + 1) q' (u) (v + 1) , 

q' n (u) = i n+1 Q' n (iu), 
1o( u ) = 2TT- 



Hence that for n — the solutions is 



,i> - 



(1 + a)e c °' 



+ (1 - a)e- c « 



— C() cot -1 u ' 



A = 

u* + v z 
A = kbc (v + 1) 



1 



(59) 
(60) 

(61) 
(62) 

(63a) 

(63b) 
(63c) 



This solution is the magnetized version of the Bonnor-Sackfield vacuum solution or zero order Morgan-Morgan 
solution [42|, [88| . The solution is regular on the axis of symmetry but is not asymptotically flat and represents the 
exterior gravitational field of finite disklike source in the presence of a magnetic field spread throughout the space 
(infinite) . 

For the terms n = and n = 2 of series (1551) we have 



(64) 



V3 = -c cot _1 (-it) - i C2 (3v 2 - l) [(3u 2 + l)cot _1 (u) - 3tt] 



For the second term of the summation we have 

C2 



.4 



kb ^l ( u a + i) q '^ u ) [ v p 2 ( v ) _ Pl ( v )} 
5 

l -kbc 2 (u 2 + 1) v (1 - v 2 ) q' 2 (u) 

\kbc 2 v (1 - v 2 ) [3u (u 2 + 1) cot _1 (u) - 3u 2 - 2] 



(65) 



Finally, the solutions is 



A = kbc (« + !) + -kbc 2 v (1 - v 2 ) [3it (u 2 + l) cot _1 (u) - 3u 2 - 2] , 



A = - 



16fc 2 



4p 2 (9u z v z - vf + v A - 1) cot _i (u)' 



+ -c 2 u(l - v 2 ){21c 2 u 2 v 2 - 3c 2 u 2 + 21c 2 u 2 - 5c 2 + 8c ) cot 1 (u) 
8 

l , N2i fu 2 +V 2 

+ ^ c 2 - c ) In — — j 
2 \ 1 + u A 

-^rc 2 (l - v 2 )(12c 2 v 2 + 27c 2 u 2 v 2 - 3c 2 u 2 - 4c 2 + 16c ), 
16 



(66a) 



(66b) 



and again ip is given by (|15a|) . In the absence of magnetic field b = and cq = c 2 , this solution corresponds 
the first order Morgan-Morgan vacuum solution. Therefore the case 6^0 corresponds to a magnetized first order 
Morgan-Morgan solution. 
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D. Solutions in generalized spheroidal coordinates 



We can reunite the three general solutions presented above in a single family by considering generalized spheroidal 
coordinates (£,/?) which are a generalization of the spherical coordinates, prolate and oblate. The asymptotically fiat 
general solution of Laplace's equation (|14l) can be written in this coordinates 



(67) 



n=0 



where c„ are constants and e = The expressions corresponding to the metric functions and the magnetic potential 
are obtained in the same way as the previous cases, that is 



A{e,rj) = / A, v dT) 



-l 



So the magnetic potential A is 

oo 

a = -kb (e a 2 ) £ J —^—^Q' n{£ ) [ V p n ( V ) - p n ^m, 

and the metric function A is given by [8(J l89j 

oo 



n,m— 



where Q^(£) is the total derivate of Q n (£,) with respect to £ and 



T m ™ = - In 
2 



1 



e 2 — rj 2 



(k n + k m 2k n k m ^ In 



S + 7] 



e- 1 



+ (e 2 - 1) 



(1 ^n) <Sm + k n S m +i 



m + 1 
1 



(P m -(-i) m QL) 



n + 1 



Am : nQmQr. 



with 



1 for I even 
for I odd 



and where A n m , B nm , C n m and S n are 



An,m ~t~ A m n P n Pj 



-1) 



f 

J-x 

j nP' m Pndv, 



n,m — ^n.m — 
(n-2)/2 

fc=0 



1 



1 



n - 2k n-2k-l 



P, 



n-2k-l 



+ (-l) n+L )Q'n-2 k -l 



(68a) 
(68b) 

(69) 
(70) 



(71) 



(72) 

(73a) 
(73b) 
(73c) 

(73d) 

Finally, the function tp is given by equation (|15al) . and with 6^0 the solutions are asymptotically flat for n > 1, 
i.e., if the coefficient cq is chosen equal to zero. So in presence of magnetic field all discussed solutions are either 
singular, not asymptotically flat or massless. 
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III. RELATIVISTIC THIN DISKS 



In order to obtain a solution of the Einstein-Maxwell equations (|2a|) - (f2q) representing a thin disk at z = with 
current, we assume that the components of the metric tensor and of the electromagnetic potential are continuous 
across the disk, but with first derivatives discontinuous on the plane z = 0, with discontinuity functions 

bab =9ab,z\z =0 + ~ 9ab,z\ s=0 - = 2 9ab,z\^ =0+ , (74) 

a a =A a>z \ z=Q+ - A a>z \^ o _ = 2A a , z \ z=Q+ . (75) 



By using the distributional approach [90 92j or the junction conditions on the extrinsic curvature of thin shells [931.19 
the Einstein-Maxwell equations yield an energy-momentum tensor T a b = T^ + T^ = T° b n + Q a b S(z), and a current 
density J a = j a S(z) = —e 2 ^~ A ^a a S(z), where S(z) is the usual Dirac function with support on the disk, T^ m is the 
electromagnetic tensor defined in Eq. (|2bj) . j a is the current density on the plane z = 0, and 

Qb = l{b az S z b ~ b z *6% + g az ¥ b - g»l$ + b c c (g zz S a b - g az S z b )} 

is the distributional energy-momentum tensor. The "true" surface energy-momentum tensor (SEMT) of the disk, S a b, 
and the "true" surface current density, j a , can be obtained through the relations 



Sab = J T^ at ds n = e A ^ Q ab , (76a) 

ja = J Ja d Sn = e A ~^ Ja , (76b) 

where ds n — y/g Z z dz is the "physical measure" of length in the direction normal to the disk. For the metric ([T]), the 
nonzero components of <S* are 

5 °= 2f>- A {A, z - 2</>, z }, (77a) 

Sl= 2e^- A A, z , (77b) 
and the only nonzero component of the current density j a in the magnetostatic case is 

h = -h e *~ AA > z > (78) 



where all the quantities are evaluated at z = + . 

In order to give physical significance to the components of the energy-momentum tensor and the electric current 

, t 

3(a) 

V a =e-+ (1,0,0,0), (79a) 



density j a we project them onto the orthonormal tetrad e( Q ) & = {V b , W b , X°, Y b }, where 



W a = — (0,1,0,0), (79b) 
P 

X a = e^- A (0, 0,1,0), (79c) 

Y a =e lJ >- A (0, 0,0,1). (79d) 

In terms of this tetrad (or observer with four-velocity V a ) the surface energy density e, the azimuthal pressure p v , 
and the azimuthal current density j of the disk are given by 

e = S°, Pv = Si, j = W 1 j 1 . (80) 

Finite thin disks can be obtained introducing oblate spheroidal coordinates, which are naturally adapted to a disk 
source. These solutions, in the vacuum and static case, correspond to the Morgan and Morgan solutions [25j . In the 
case of the magnetized Morgan-Morgan solutions, following the reference (25^ . for n even the metric functions are 
continuous across the disk but its first derivatives are discontinuous in the direction normal to the disk, which can be 
interpreted as a finite thin disk located at z — and < p < 1. However, since for n even the function A is a odd 
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polynomial of v, the opposite occurs with the magnetic potential, i.e, it is discontinuous across the disk but its first 
derivatives are continuous in the direction normal to the disk, so that electric current density also is zero on the disk 
and in consequence the source of the magnetic field is non planar but of a different origin such as a remnants or fossil 
magnetic field [62], or can come from external sources, such as the presence of a nearby magnetars or neutron stars. 
Indeed on the disk the Maxwell equations dbF ab — 4irJ a , where 'bar' denotes multiplication by y/—g, are given by 

- A-Kj v S(z) = d z g zz A, z + d p g pp A^ (81) 

where S(z) is the usual Dirac Function with support on the disk. Integrating through the disk we obtain 

rz=0+ r z=0 



±k% = f + d z g zz A z dz + f + d p g pp A <p dz 

3+ 



I 2 = 



9 ZZ Az\lZ + +d p g"Pd p l Z ° + Adz 
= 0, (82) 

where the first term on the right-hand side vanishes due to the continuity of the metric and A >z and the second term 
from discontinuity of A, or in other words as A is an odd polynomial function of v its integral through the disk is 
even and hence continuous, so that to the evaluate at the limits of integration expression vanishes. Thus we can 
interpret this solutions as the exterior gravitational field of a finite dislike source immersed in a magnetic field. For 
completeness, we analysis the electrostatic case, i.e., for $ real. The electric potential <fi is given by 

P[ ' (83) 



(1 + a)e~^ + (1 - a)ev 



where a = y/l + p 2 , being p the parameter that controls the electric field, and the only non-zero component of the 
current density j a on the plane z = is 

.it = -^~ A ^ (84) 
whereas in terms of the tetrad (|79p the electric charge density a is given by 

a = -V°j . (85) 

We see that just as metric the electric potential <f> is a even function of v and in consequence the current density in 
the disk is nonzero. Therefore, this solutions can be interpreted as the gravitational field of a finite charged disk. 

In order to study the behavior of the main physical quantities associated with these disks we perform a graphical 
analysis of them for magnetized first order Morgan-Morgan finite disks. In Fig. [I] we show the energy density e and 
the azimuthal pressure p v with c = c 2 = 0.4 and for values of magnetic field parameter b — (dashed curves), 0.5, 
1, and 2 (bottom curves), as functions of p. We see that the energy density presents a maximum at p = and then 
decreases with p. We also see that the presence of magnetic field decreases the energy density at the central region 
of the disk and later increases it. We can observe that the pressure increases rapidly as one moves away from the 
disk center, reaches a maximum and later decreases. We also observe that the magnetic field decreases the pressure 
everywhere on the disk. The graph also show that the disk's surface energy density is always positive in concordance 
with the weak energy condition, as well as the stress in azimuthal direction (pressure). The strong energy condition, 
e+Ptp > 0, is also satisfied. These properties characterize a fluid made of matter with the usual gravitational attractive 
property. In addition, in Fig. QTc) we show, in the electrostatic case, the surface electric charge density a for values of 
electric field parameter p = 0.2 (dashed curve), 0.5, and 1 (top curve) and the same values of cq and c-i. We observe 
a greater concentration of electric charge near the rim of the disk and that this increases as we increase the electric 
field. 

Exact solutions which represent the field of a disk also can be obtained using the well known "displace, cut and 
reflect" method that was first used by Kuzmin |95l| and Toomre [96( to constructed Newtonian models of disks, and 
later extended to general relativity [33|, H3, [H, |37I |. Given a solution of the Einstein-Maxwell equation, this procedure 
is mathematically equivalent to apply the transformation z — > \z\ + z$, with zq constant, on that solution. However, 
this disks are essentially of infinite extension and the field not correspond exactly to the metric with which we started. 
These solutions can be called in our case magnetized Weyl type solutions. The reflection symmetry with respect to 
the plane of disk implies that in this case A is continuous across the disk and its first derivative discontinuous in the 



14 



direction normal to the disk, so that electric current density is nonzero on the disk and in consequence the source of 
the magnetic held is planar. 

In Figs. [2] - 14] we show, as functions of p, the energy density e, the azimuthal pressure p v and surface azimuthal 
electric current density j for magnetized Weyl type disks of infinite extension corresponding to the lineal combination 
of the first and second terms of the series in the three families of solutions considered. Since the surface energy density 
decreases rapidly one can to dehne a cut off radius and, in principle, to consider these disks as finite. Anyway all the 
physical quantities present a similar behavior to the previous case. 



IV. MOTION OF CHARGED PARTICLES AROUND MAGNETIZED WEYL FIELDS 



The relativistic Lagrangian for a test particle in presence of a gravitational and magnetic field is given by 

C = l -9abx a x h + eA<p, (86) 

where e = e/m is the specific charge of the particle and the overdot denotes derivate with respect to the proper time 
s. For magnetized Weyl fields we have two constants of motion 

E = -p t /m = -g tt i, (87) 
L = p v /m = g vv <p + eA, (88) 

where E represents the relativistic specific energy and L the specific angular momentum. 

For circular orbits and from symmetry of the field the equation for the electrogeodesic motion of the particle is 
given by 

\gab, P u a u b = -eF pa u a (89) 

where u a is the 4-velocity of particles with respect a the coordinates frame. For equatorial circular orbits u a has 
components u a = w°(l,w,0, 0) where lu — u 1 /u° — -£ is the angular velocity of the test particles, and in the case 
magnetostatic the equation of motion reads 



^M°(5 w ,pW 2 + gtt, P ) = -eA !P u>, (90) 
where u° obtains normalizing u a , that is requiring g a i,u a u b = — 1. Thus with 

(u°f = J— - (91) 

g w u z + g tt 

the equation of the electrogeodesic takes the form 

A(uj 2 ) 2 + Buj 2 + C = 0, (92) 

where 

A = 9 2 W , P + 4e 2 g vv A 2 p 

= 4p 2 e- 2 ^ [(1 - piPJ 2 e- 2 ^ + e 2 A 2 p ] , (93a) 
B = 2g tt . p g v>VtP + Ae 2 g tt A 2 p 

= -4[2p^ P (l - + e 2 A 2 p e 2 ^ , (93b) 
C = gl p = ^ 2 p e^. (93c) 
Therefore, the angular velocity uj is given by 

2 -B±y/D 

where 

D = B 2 - 4AC = \<ol 2 A 2 e 2 ^\\pi), p (\ - 2pV, P ) + e 2 A 2 e 2 ^] > 0, (95) 
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and positive sign corresponds to the direct orbits (or co-rotating with L > 0) and the negative sign to the retrograde 
orbits (or counter-rotating with L < 0). From (pTTj) and we find 

E = -g tt u°, (96) 
L = g vv uju + eA. (97) 

With respect to the orthonormal tetrad (|79|) the 3- velocity has components 

»W = 1»!L (98) 

For equatorial circular orbits the only nonvanishing velocity components is given by 

(^) 2 ^ c 2 = -^ W 2 = p 2 e- 4 ^ 2 , (99) 
9tt 

which is the circular velocity of the particle as seen by an observer at infinity. In fact when p — > oo, tp — > and 
v 2 = p 2 ui 2 = pjf-, where <fr is the Newtonian gravitational potential, and v c is the circular velocity which represents 
the velocity of a test particle in a circular orbit at radius p. For neutral particles e = we have 

< = (100) 



with 



,l + o)e-^-(l-o)e^ 
(l + o)e-^ + (1 - a)e^ 



(101) 



In terms of the circular velocity v c the specific angular momentum can be written as 

To analyze the stability of circular orbits on the equatorial plane we use an extension of Rayleigh criteria of stability 
of a fluid at rest in a gravitational field |68i470j . The motion equation (|90p can be cast as a balance equation 



9 pp 9u, P E 2 Q gPPg^ L-eA) 

2 9 l +£9 ^ 2 9{pv 



2 + e 5 "M, pWU ° = - 2 ™^ ^ (103) 



where the term first on the left-hand side represents the gravitational force F g , the term second the the Lorentz force 
Fl, and the term on the right-hand side the centrifugal force F c (p) = F(p, L(p)) acting on the particle. So we have a 
balance between the total force F(p) = F g + Fl and the centrifugal force. We now consider the particle to be initially 
in a circular orbit with radius p = po and we slightly displace it to a higher orbit p > pq. The angular momentum 
of particle remains equal to its initial value Lq — L(po) which implies that the centrifugal force in its new position is 
F c (p,Lq). In order that the particle returns to it initial position must be met that F(p) > F c (p, Lq), but according 
to the balance equation (|103[) F(p) = F c (p, L) so that F c (p, L) > F c (p, Lq), and hence (L — eA) 2 > (L — eA) 2 . Using 
the expresion for L (1UT1) and defining la function h — g^^uiu , follows that h(p) 2 > h(p ) 2 . Note that the quantity h 
has the same form that the specific angular momentum in the vacuum. By doing a Taylor expansion of h 2 (p) around 
P = Po one finds that the condition of stability for a circular orbit is 

hh tP > 0, (104) 

or, in other words, 

h 2 p > 0. (105) 

In order to study the behavior of these physical quantities for the previous solutions again we perform a graphical 
analysis of them for the zeroth and second terms of the series in the three families of solutions considered. In all cases 
we take e = 1. We first analyze the magnetized Weyl solutions in spherical coordinates. In this case A p = 0, so that 
the velocity is given by expression (jlOOl) . In Figure [S] we show the circular velocity curves v 2 and h 2 for test particles 
with Co = 0.5 and C2 = 1 (figures [HJa) and , and cq = and c-i = —0.5 
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(figures[5jc) andHfc?)), for values of magnetic field parameter 6 = (dashed curves), 1, 2, 3 (top curves), as functions 
of p. In first case we see that the tangential velocity increase initially from certain p — p , reaches a maximum and 
then falls to zero, and is always a quantity less than the velocity of light. We also observer that inclusion of magnetic 
field make these orbits less relativistic. Meantime we see that the quantity h 2 always is a increasing function of p 
that corresponds to stable circular orbits. In the second case we find that v 2 increases from zero at infinity to 1 (the 
velocity of light) at circular photon orbit p = p p h, and that h 2 always is a decreasing monotonous function of p which 
means that this orbits are instable against radial perturbation. 

We now analyze the magnetized Erez- Rosen solutions. In this case also A tP = 0. So in figure [5] we show the circular 
velocity curves v 2 and h 2 for test particles with Co = 0.5 and c 2 = 1 (figures [B^a) and|U[6)), and Co = and c 2 = —0.5 
(figures |6jc) andlUJd)), for values of magnetic field parameter 6 = (dashed curves), 0.5, 1, 2 (top curves), as functions 
of p. In first case we see that the tangential velocity increase initially, reaches a maximum and then falls to zero, 
being always a quantity less than the velocity of light. We also observer that inclusion of magnetic field make these 
orbits less relativistic. The quantity h 2 always is a increasing function of p that corresponds to stable circular orbits. 
In the second case we find that v 2 also increase initially but not from zero and that only the central regions are stable. 

In Figures [7J and [5] we show for magnetized first order Morgan-Morgan fields the circular velocity curves v 2 and h 2 
for test particles with cq — 0.5, c 2 = 1 (figures [7]) and with Co = 0, c 2 = — 1 (figures [5]) for values of magnetic field 
parameter 6 = 0,1, and 1.4. We find inside of the disk A :P ^ so that in this case we have two values for the angular 
velocity w±, which corresponds to the direct and retrograde motions of particles. Figures on the left side correspond 
to the inside of the disk and to the direct motion of particles. For the retrograde motion of the particles we obtain 
a similar behavior. For all values of parameters we find regions where circular orbits are possible, that is, regions 
where the velocity of particles is always a quantity less than the velocity of light. We also observer that inclusion of 
magnetic field can makes these orbits less relativistic. We find strong change in the slope of h 2 at certain values of p 
which means that there is a strong instability there, also regions with negative slope where the orbits are unstable, 
but always we find regions where h 2 ^ is a increasing monotonous function of p that corresponds to stable circular 
orbits. We also observe that the increase of magnetic field can make stable these orbits against radial perturbation. 

In order to compare the behavior of the se q uantities with other known magnetized solutions we now analyze a Kerr 
type solution or magnetic dipole solution [74| . In this case also we have two values for the angular velocity, ui±. For 
the direct motion of the particles, in figure [9] we present EJa) the circular velocity curves v 2 + for 6 = (dashed curve), 
0.5 (curve with points and lines), 1.5, 2, 2.5, and 3 (bottom curve), and h? + for[Sj6) 6 = (dashed curve), 0.5, Etc) 
6 = 1.5 (dashed curve), 2, 2.5, and 3 (bottom curve), as functions of p. We see that in absence of magnetic field 6 = 
the circular velocity increases from zero at infinity to 1 (the velocity of light) at circular photon orbit p = p p h, but as 
we increase the magnetic field it reaches a maximum and then falls to zero, being always less than the velocity of light. 
We also see that inclusion of magnetic field make these orbits less relativistic. Regarding the function h 2 we find for 
6 = and 6 = 0.5 (or weak fields) strong change in the slope of h 2 + at certain values of p which means that there is a 
strong instability there, also regions with negative slope where the orbits are unstable, but after certain values of p we 
find that h 2 ^ is a increasing monotonous function of p that corresponds to stable circular orbits. We also observe that 
the increase of magnetic field can make stable these orbits against radial perturbation. For the retrograde motion of 
the particles we obtain a similar behavior. 



V. DISCUSSION 

Using the well-known Ernst's method were constructed three families of general magnetostatic axisymmetric exact 
solutions of Einstein-Maxwell equations in spherical coordinates, prolate, oblates, and also in generalized spheroidal 
coordinates which is a generalization of the previous systems. Except the zero order, the constructed solutions are 
asymptotically flat but with total mass zero. Hence, all discussed solutions with magnetic field are either singular, not 
asymptotically flat or massless. In the three cases considered we show explicitly the solutions corresponding to the 
zeroth and second terms of the series of the solutions, among them a magnetized Erez- Rosen metric and a magnetized 
Morgan-Morgan metric, which was interpreted as the exterior gravitational field of a finite dislike source immersed 
in a magnetic field, and we analyze the material properties of the disk such as the surface energy density and the 
azimuthal pressure. From them we also study relativistic models of thin disks of infinite extension. We find always 
values of parameters for which all the energy conditions are satisfied. 

We also analyzed the electrogeodesic equatorial circular motion of test particles for the sum of the first and second 
terms of series of solutions and we also discuss their stability. For all values of parameters we find regions where 
circular orbits are possible, that is, regions where the circular velocity of particles is always a quantity less than the 
velocity of light. We also observer that inclusion of magnetic field can makes these orbits less relativistic. We find 
strong change in the slope of the function h 2 at certain values of the radial coordinate p which means that there is 
a strong instability there, also regions with negative slope where the circular orbits are unstable, but always we find 
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values of parameter for which h 2 is a increasing monotonous function of p that corresponds to stable orbits. We also 
observe that the increase of magnetic field can make stable these orbits against radial perturbation. For magnetized 
Morgan-Morgan fields was found that inside of disk the presence of magnetic field provides the possibility of to find 
relativist charged particles moving in both direct and retrograde direction. 

Finally in order to compare the behavior of these physical quantities with other known magnetized solutions, we 
also analyze the electrogeodesic motion of test particles and their stability for a Kerr-type solution (magnetic dipolc 
solution). We found for some value of the parameter a similar behavior to the previous cases. 



[1] A.P. Lobanov, A&A 330, 79 (1998). 

[2] F.A. Aharonian, MNRAS 332, 215 (2002). 

[3] A.F. Zakharov et al, MNRAS 342, 1325 (2003). 

[4] J.S. Greaves et al, Nature 404, 732 (2000). 

[5] Y.T. Liu, S.L. Shapiro and B.C. Stephens, Phys. Rev. D 76, 084017 (2007). 
[6] S. Akiyama et al, ApJ 584, 954 (2003). 

[7] W.H.T. Vlemmings, RJ. Diamond and H. Imai, Nature 440, 58 (2006). 

[8] D.L. Meier, S. Koide and Y. Uchida, Science 291, 84 (2001). 

[9] G.V. Ustyugova et al, ApJ 516, 221 (1999). 
[10] R. Krasnopolsky, Z.Y. Li and R.D. Blandford, ApJ 595, 631 (2003). 
[11] D.N. Aguilera, J.A. Pons and J.A. Miralles, ApJ 673, L167 (2008). 
[12] P. Hennebelle and S. Fromang, A&A 477, 9 (2008). 
[13] Z. Medin and D. Lai, MNRAS 382, 1833 (2007). 
[14] X.D. Li and E.P.J, van den Heuvel, ApJ 513, L45 (1999). 
[15] G.G. Pavlov and V.G. Bezchastnov, ApJ 635, L61 (2005). 
[16] A.I. Ibrahim, J.H. Swank and W. Parke, ApJ 584, L17 (2003). 

[17] K. Schwarzschild, Sitzungsberichte der Kniglich Preussis- chen Akademie der Wissenschaften 1, 189 (1916). 

[18] K. Schwarzschild, Sitzungsberichte der Kniglich Preussis- chen Akademie der Wissenschaften 1, 424 (1916). 

[19] R.P. Kerr, Phys. Rev. Lett. 11, 237 (1963). 

[20] E.T. Newman et al, J. Math. Phys. 6, 918 (1965). 

[21] H. Reissner, Ann. Phys. Berlin 50, 106 (1916). 

[22] G. Nordstr m, Verhandl. Koninkl. Ned. Akad. Weten-oschap. 20, 1231 (1918). 

[23] J. Polanco, P. Letelier, and M. Ujevic , Phys. Rev. D 78, 024026 (2008). 

[24] W. A. Bonnor and A. Sackfield, Commun. Math. Phys. 8, 338 (1968). 

[25] T. Morgan and L. Morgan, Phys. Rev. 183, 1097 (1969). 

[26] L. Morgan and T. Morgan, Phys. Rev. D 2, 2756 (1970). 

[27] G. A. Gonzalez and P. S. Letelier, Class. Quantum Grav. 16, 479 (1999). 

[28] D. Lynden-Bell and S. Pineault, Mon. Not. R. Astron. Soc. 185, 679 (1978). 

[29] A. Chamorro, R. Gregory, and J. M. Stewart, Proc. R. Soc. London A413, 251 (1987). 

[30] PS. Letelier and S. R. Oliveira, J. Math. Phys. 28, 165 (1987). 

[31] J. P. S. Lemos, Class. Quantum Grav. 6, 1219 (1989). 

[32] J. P. S. Lemos and P. S. Letelier, Class. Quantum Grav. 10, L75 (1993). 

[33] J. Bicak, D. Lynden-Bell, and J. Katz, Phys. Rev. D 47, 4334 (1993). 

[34] J. Bicak, D. Lynden-Bell, and C. Pichon, Mon. Not. R. Astron. Soc. 265, 126 (1993). 

[35] G.A. Gonzalez and O. A. Espitia, Phys. Rev. D 68, 104028 (2003). 

[36] J. Bicak and T. Ledvinka, Phys. Rev. Lett. 71, 1669 (1993). 

[37] G. A. Gonzalez and P. S. Letelier, Phys. Rev. D 62, 064025 (2000). 

[38] T. Ledvinka, J. Bicak, and M. Zofka, in Proceeding of 8th Marcel-Grossrnann Meeting in General Relativity, edited by T. 

Piran (World Scientific, Singapore, 1999) 
[39] G. Garcfa-Reyes and G. A. Gonzalez, Brazilian Journal of Physics 37, no. 3B, 1094 (2007). 
[40] P. S. Letelier, Phys. Rev. D 60, 104042 (1999). 

[41] J. Katz, J. Bicak, and D. Lynden-Bell, Class. Quantum Grav. 16, 4023 (1999). 

[42] G. Garcia R. and G. A. Gonzalez, Phys. Rev. D 69, 124002 (2004). 

[43] D. Vogt and P. S. Letelier, Phys. Rev. D 70, 064003 (2004). 

[44] G. Garcfa-Reyes and G. A. Gonzalez, Class. Quantum Grav. 21, 4845 (2004). 

[45] G. Garcfa-Reyes and G. A. Gonzalez, Phys. Rev. D 70, 104005 (2004). 

[46] J. L. Synge, Relativity: The General Theory. (North-Holland, Amsterdam, 1966). 

[47] G. Neugebauer and R. Meinel, Phys. Rev. Lett. 75, 3046 (1995). 

[48] C. Klein, Class. Quantum Grav. 14, 2267 (1997). 

[49] C. Klein and O. Richter, Phys. Rev. Lett. 83, 2884 (1999). 

[50] C. Klein, Phys. Rev. D 63, 064033 (2001). 

[51] J. Frauendiener and C. Klein, Phys. Rev. D 63, 084025 (2001). 



18 



C. Klein, Phys. Rev. D 65, 084029 (2002). 
C. Klein, Phys. Rev. D 68, 027501 (2003). 

C. Klein, Ann. Phys. (Leipzig) 12, 599 (2003). 

D. Vokrouhlicky and V. Karas, Astron. Astrophys. 243, 165 (1991). 
S. Sengupta, Int. J. Mod. Phys. D 6, 591 (1997). 

F. de Felice and F. Sorge, Class. Quantum Grav. 20, 469 (2003). 

J. Kovaf, O. Kopacek, V. Karas and Z StuchHk, Class. Quantum Grav. 27, 135006 (2010). 

J. Frank , A. King and D. Raine, Accretion Power in Astrophysics (Cambridge: Cambridge University Press, 1992). 

G. Belvedere, Proc. European Physical Society Study Conference, Accretion Disks and Magnetic Fields in Astrophysics 
(Dordrecht: Kluwer, 1998). 

A. N. Aliev and N. Ozdemir, Mon. Not. Roy. Astron. Soc. 336, 241 (2002). 
J. Braithwaite and H. C. Spruit, Nature 431, 819 (2004). 

V. C. Rubin, J. A. Graham and J. D. P Kenney. Ap. J. 394, L9, (1992). 

H. Rix, M. Franx, D. Fisher and G. Illingworth. Ap. J. 400, L5, (1992). 
F. Bertola et al. Ap. J. 458, L67 (1996). 

C. Struck, Phys. Rep. 321, 1 (1999). 

R. Ciri, D. Bettoni, and G. Galletta, Nature 375, 661 (1995). 
Lord Rayleigh, 1917, Proc. R. Soc. London A, 93, 148 

L. D. Landau and E.M. Lifshitz, Fluid Mechanics(Addison- Wesley, Reading, MA, 1989). 
P. S. Letelier, Phys. Rev. D 68, 104002 (2003). 
F. J. Ernst. Phys. Rev. 167, 1175 (1968). 

F. J. Ernst. Phys. Rev. 168, 1415 (1968). 

D. Kramer, H. Stephani, E. Herlt, and M. McCallum, Exact Solutions of Einsteins's Field Equations (Cambridge University 
Press, Cambridge, England, 1980). 

W. B. Bonnor, Z. Phys. 190, 444 (1966). 

Bateman, H. Higher Transcendental Functions. McGraw Hill, New York, vol. 2 
Arfken, G. Mathematical Methods for Phystcists. Academic Press. 
S. M. Scott and P. Szekeres, Gen. Rel. Grav., 18, 557 (1986). 
S. M. Scott and P. Szekeres, Gen. Rel. Grav., 18, 571 (1986). 

J. M. Bardeen, B. Carter, S. W. Hawking, Commun. Math. Phys., 31, 161 (1973). 
H. Quevedo, Phys. Rev. D 39, No. 10, 2904 (1989). 
Chazy J 1924 Bull. Soc. Math., France 52 17. 
Curzon H E J 1924 Proc. London Math. Soc. 23 477. 

S. Chandrasekhar, The Mathematical Theory of Black holes (Oxford Publications, New York, 1992), p. 623. 

D. M. Zipoy. J. Math. Phys. 7, 1137 (1966). 

B. H. Voorhees. Phys. Rev. D 2, 2119 (1970). 

K. Schwarzschild. Sitzungsber. Preuss. Akad. Wiss., 189 (1916). 

G. Erez and N. Rosen, Bull. Res. Counc. Isr. 8F, 47 (1959). 

W. A. Bonnor and A. Sackfield, Comm. Math. Phys. 8, 338 (1968). 
G. A. Gonzalez and J. F. Ramos, Rev. Col. Ffs. 33, No. 2, 118 (2001). 
A. Papapetrou and A. Hamouni, Ann. Inst. Henri Poincare 9, 179 (1968) 
A. Lichnerowicz, C.R. Acad. Sci. 273, 528 (1971) 
A. H. Taub, J. Math. Phys. 21, 1423 (1980) 

E. Israel, Nuovo Cimento 44B, 1 (1966) 
E. Israel, Nuovo Cimento 48B, 463 (1967) 

G. G. Kuzmin 1956, Astron. Zh., 33, 27 (1956) 
A. Toomrc, Ap. J., 138, 385 (1962) 



19 




P P P 

(a) (b) (c) 

FIG. 1. (a) The surface energy density e and (b) the azimuthal pressure p v for magnetized first order Morgan-Morgan finite 
disks with Co = C2 = 0.4 and for values of magnetic field parameter 6 = (dashed curves), 0.5, 1, and 2 (bottom curves), as 
functions of p. (c) The surface electric charge density a (electrostatic case) for p — 0.2 (dashed curve), 0.5, and 1 (top curve) 
and the same values of Co and c%. 
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(a) (b) (c) 

FIG. 2. (a) The surface energy density e and (b) the azimuthal pressure p v for magnetized weyl type infinite disks in spherical 
coordinates and second order with zq = 1, eo = 1, C2 = 0.5 and for values of magnetic field parameter 6 = (dashed curves), 
0.5, 1, and 2 (bottom curves), as functions of p. (c) The azimuthal current density j for b = 0.5 (dashed curve), 1, and 2 (top 
curve) and the same values of zo , Co and ci . 
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FIG. 3. (a) The surface energy density e and (b) the azimuthal pressure p v for magnetized Erez-Rosen type infinite disks with 
zo — 2, Co = C2 = 1 and for values of magnetic field parameter 6 = (dashed curves), 1, 2, and 3 (bottom curves), as functions 
of p. (c) The azimuthal current density j for 6=1 (dashed curve), 2, and 3 (top curve) and the same values of zq, cq and C2- 




FIG. 4. (a) The surface energy density e and (b) the azimuthal pressure for magnetized first order Morgan-Morgan type 
infinite disks with zo = Co = C2 = 1 and for values of magnetic field parameter b = (dashed curves), 0.5, 1, and 2 (bottom 
curves), as functions of p. (c) The azimuthal current density j for 6 = 0.5 (dashed curve), 1, and 2 (top curve) and the same 
values of zo, Co and C2- 
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FIG. 5. For magnetized weyl fields in spherical coordinates and second order we plot, as functions of p, the circular velocity v 2 
and h 2 for test particles with Co = 0.5 and C2 = 1 (figures (a) y (b) ), and co = and C2 = — 1 (figures (c) y (d) ), for values of 
magnetic field parameter b = (dashed curves), 1, 2, 3 (top curves). 
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FIG. 6. For magnetized Erez-Rosen fields we plot, as functions of p, the circular velocity v% and h 2 for test particles with 
co = 0.5 and ca = 1 (figures (a) y (6) ), and Co = and C2 = —0.5 (figures (c) y (d) ), for values of magnetic field parameter 
6 = (dashed curves), 0.5, 1, 2 (top curves). 
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FIG. 7. For magnetized first order Morgan-Morgan fields we plot, as functions of p, the circular velocity v 2 and h 2 for test 
particles with Co = 0.5 and C2 = 1, for values of magnetic field parameter 6 = (dashed curves), 0.5, 1, and 1.4 (top curves). 
Figures on the left side correspond to the inside of the disk and to the motion direct of particles. 
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FIG. 8. Also for magnetized first order Morgan-Morgan fields we plot, as functions of p, the circular velocity v\ and h 2 for test 
particles with Co = and C2 = —1, for values of magnetic field parameter 6 = 0, 0.5, 1, 1.4 (dashed curves). Again figures on 
the left side correspond to the interior region of the disk and the motion direct of particles. 
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FIG. 9. For a Kerr-type solution (magnetic dipole solution) we plot, as functions of p, (a) the circular velocity vf+ for test 
particles with 6 = (dashed curve), 0.5 (curve with points and lines), 1.5, 2, 2.5, and 3 (bottom curve), and h\ for (b) b = 
(dashed curve), 0.5, (c) b = 1.5 (dashed curve), 2, 2.5, and 3 (bottom curve). 



